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ABSTRACT 

The recent discovery, by the Chandra satellite, that jets of blazars are strong X-ray 
emitters at large scales (0.1-1 Mpc) lends support to the hypothesis that emitting 
plasma is still moving at highly relativistic speeds on these scales. In this case in 
fact the emission via inverse Compton scattering off cosmic background photons is 
enhanced and the resulting predicted X-ray spectrum well accounts for the otherwise 
puzzling observations. Here we point out another reason to favour relativistic large 
scale jets, based on a minimum power argument: by estimating the Poynting flux and 
bulk kinetic powers corresponding to, at least, the relativistic particles and magnetic 
held responsible for the emission, one can derive the value of the bulk Lorentz factor 
for which the total power is minimized. It is found that both the inner and extended 
parts of the jet of PKS 0637-752 satisfy such a condition. 
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1 INTRODUCTION 

Most of the radiation observed from blazars (BL Lac objects 
and flat spectrum radio quasars) originates in relativistic 
jets, in which the emitting plasma flows with a bulk Lorentz 
factor F ~ 10 at small viewing angles 6. Since the observed 
radiation is beamed in the direction of motion the large ap- 
parent luminosities are much reduced in the comoving frame, 
by a factor T 4 (if 8 ~ 1/r) in the case of a moving blob. On 
the other hand, the most important power content of the 
jet is not in the form of internal random particle energy 
giving rise to radiation, but in the form of electromagnetic 
and kinetic powers associated with the bulk motion. This 
is indicated also by the power budget of the extended radio 
structures and radio lobes, which require an average supply 
which is greater than the luminosity generated by the jet 
at all scales. However, the jet does not need to remain rela- 
tivistic for all its length, and indeed it was believed that it 
slows down, without losing much of its bulk kinetic energy 
content. In fact the detection of large-scale, symmetric jets 
has been considered as an indication that on such scales the 
plasma is moving at most with mildly relativistic velocities 
in low power FR 1 sources (Fanaroff & Riley 1974), while 
the jet to counterjet flux ratio requires /3c ~0.6-0.8 in FR II 
large-scale jets (e.g. Laing 1993; Bridle et al. 1994; Wardle 
& Aaron 1997; Hardcastle et al. 1999). However, as these 
authors (and Laing et al. 1999) pointed out, there is the pos- 
sibility of a velocity structure across the jet, whose plasma 
could move highly relativistically close to the axis ( "spine" ) 



and mildly relativistic at the border ( "slow layer" ) . As found 
by Chiaberge et al. (2000), Celotti et al. (2001), such a ve- 
locity structure can also help to unify radio-galaxies seen at 
large viewing angles and blazars. 

Observations of large scale jets emitting copious 
amounts of X-rays also suggest that jets might maintain 
ultra-relativistic speed up to the largest scales. We consider 
as a specific case that of the source PKS 0637-752, observed 
by Chandra, for which observational information is available 
on both large and small scales. Chandra observations reveal 
relatively intense X-ray emission quite far from the nucleus 
(Chartas et al. 2000). The spatially resolved X-ray images 
show knots of emission similar to those seen in radio maps 
(Tingay et al. 1998) and optical images (HST; Schwartz et 
al. 2000). VSOP observations detected superluminal motion 
in the jet associated with the source (Lovell 2000), setting a 
lower limit on the (inner, i.e. pc scale) bulk Lorentz factor 
r > 17.5 and an upper limit on the viewing angle 9 < 6.4 
degrees (assuming Ho = 50 km s^ 1 Mpc -1 and qo = 0). Fur- 
thermore the small scale jet appears to be well aligned with 
the kpc scale one: by de-projecting the X-ray map using an 
angle of 6 degrees the observed X-ray emission extends up 
to ~1 Mpc from the nucleus. 

As discussed by Chartas et al. and Schwartz et al., on 
large scales the radio spectrum can be explained as syn- 
chrotron emission, as also supported by the detection of 
strong linear radio polarization, while the X-ray emission 
exceeds by a large factor (~ 2 orders of magnitude) the 
extrapolation from the radio-optical spectrum. If the X- 
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ray flux is synchrotron self-Compton (SSC) radiation from 
plasma moving at most at mildly relativistic speed, then 
large deviations from equipartition between magnetic field 
and particle energy density and/or large inhomogeneities 
have to be invoked in the emitting region. On the other hand 
thermal emission by e.g. shocked plasma also appears to be 
unlikely, as it would imply high particle density and thus 
unobserved large rotation measures. On the other hand the 
X-ray flux can be satisfactorily interpreted as scattered Cos- 
mic Microwave Background radiation (CMB) by relativistic 
plasma moving with a bulk Lorentz factor Y ~10-15 (Ghis- 
ellini 2000; Celotti, Ghisellini & Chiaberge 2001; Tavecchio 
et al. 2000). These observations then suggest that the jet, 
or some part of it, can indeed maintain highly (bulk) rela- 
tivistic speeds from milli- to megaparsec. 

Here we introduce another argument in support of large 
bulk Lorentz factors along the entire jet extension. This is 
based on minimizing the power required to produce the ra- 
diation we observe, once we consider not only the dissipated 
radiative luminosity, but also the Poynting flux and the par- 
ticle bulk kinetic power of the jet. Note that, as the power 
in bulk motion dominates over the internal particle energy, 
the minimum power criterion does not coincide with the 
equipartition argument usually applied to radio lobes (Bur- 
bidge 1959). We show that the behaviourof the total power 
as function of the bulk Lorentz factor depends on the dom- 
inant emission process (chosen to reproduce the intrinsic 
quantities) and has a minimum for a certain value of T. 

In the following we first explicitly derive the condition 
for the minimum total power requirement (Section 2). This is 
applied to the X-ray extended emitting knot and small scale 
jet of PKS 0637-752 in Section 3. The implications for the 
energy transport in jets and some observational predictions 
are discussed in Section 4. 



2 THE MINIMUM TOTAL POWER 

The spectral energy distribution (SED) of blazars and more 
generally the emission by non-thermal plasma in relativistic 
jets can usually be represented by one or two broad com- 
ponents. These are commonly interpreted as synchrotron 
(the low energy one) and possibly inverse Compton emis- 
sion (the high energy one, if present) from a population 
of relativistic leptons. The seed photon population for the 
inverse-Compton scattering may be the synchrotron pho- 
tons themselves (synchrotron self-Compton, SSC) or pho- 
tons produced externally to the emitting plasma. 

Let us assume the simplest scenario and approximate 
the emitting region as a spherical 'blob' of plasma of size R, 
in bulk motion with Lorentz factor Y at an angle d with re- 
spect to the line of sight. This produces an intrinsic power L' a 
in synchrotron radiation, and L' c in inverse Compton radia- 
tion. These are the luminosities as measured in the comoving 
frame. The corresponding dissipated luminosities (i.e. the 
power received by a detector completely surrounding the 
source and which sees it moving with F) are enhanced by 
the factor ~ (4/3)T 2 . Instead the 'observed' luminosities L B 
and L c , estimated by multiplying the received flux by 4-7rd 2 
(cLl is the luminosity distance), are a factor 8 4 greater then 
the intrinsic ones. Here 5 = [r(l — /3cos#)] -1 is the Doppler 
factor and (5c the blob bulk velocity. 



The plasma comprises leptons (protons) of comoving 
density n (n p ), embedded in a magnetic field of strength B 
perpendicular to the direction of motion. The power carried 
as Poynting flux and the kinetic powers corresponding to the 
bulk motion of leptons and protons are respectively given by: 



-R 2 Y 2 j3cB 2 . 
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Here (7) is the average random Lorentz factor of the rela- 
tivistic leptons, measured in the comoving frame, and m p , 
m are the proton and electron rest masses, respectively. 
Protons are assumed to be cold in the comoving frame. 

We can estimate these quantities from the observed ra- 
diation. A lower limit on n e can be inferred from the syn- 
chrotron emission L s of the plasma. In fact the number den- 
sity of leptons producing the observed radiation is: 



(2) 
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where (j 2 ) is averaged over the relativistic electron distri- 
bution and or is the Thomson scattering cross section. The 
particle bulk kinetic powers can then be expressed as 
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This implies that for a given size of the emitting region, 
T, 9, observed synchrotron luminosity and spectrum (which 
determines (7) and (7 2 )), the total power Lb+L c + L p must 
have a minimum for some value of B. 

A similar approach has been adopted for the case of 
the galactic superluminal source GRS 1915+105, for which 
r and 6 are known, by Gliozzi, Bodo & Ghisellini (1999). 
In the case of blazars the viewing angle is constrained by 
selection effects to be close to 6 = 1/r, since the flux is 
strongly dimmed for larger viewing angles, while the solid 
angle becomes too small for 6 < 1/F. Under this further 
condition, the total power depends only on the product YB. 
In fact in this case 8 ~ F, and then Lb and L c + L p behave in 
opposite ways with respect to YB. The value corresponding 
to the minimum total power can be easily determined by 
setting d/d(BY)(L c + L p + L B ) = 0, to derive (BF) min : 



(BY), 
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The behaviours of the different powers and the resulting 
minimum of the total one are illustrated, as an example, in 
Fig. 1. Here we represent their values as a function of YB 
in the case of "normal" electron-proton (e-p) plasma and 
pure electron-positron (e ) pair plasma. 

Interestingly a second observational constraint allows 
us to decouple B and Y. In fact the ratio of the luminos- 
ity in the inverse Compton component to that in the syn- 
chrotron component (hereafter the Compton dominance) is 
proportional to the ratio of the relevant energy densities, in 
radiation and magnetic field. 

Consider first the case of inverse Compton on photons 



© 0000 RAS, MNRAS 000, 000-000 



RELATIVISTIC LARGE SCALE JETS AND MINIMUM POWER REQUIREMENTS 3 



48 



46 




0.0001 



0.001 



0.0 1 



TB 



Figure 1. The kinetic (in protons, L p , and in leptons, L e ) and 
magnetic powers as functions of FB, the product of the the bulk 
Lorentz factor and the magnetic field. The viewing angle is as- 
sumed such that 5 = r. The axis on the top refers to values of T 
for a representative value of B. 



produced externally to the jet. The ratio between the Comp- 
ton and synchrotron luminosities can be expressed as 

L t = (5) 

where the radiation energy density U^t (~ r 2 ?7 ex t for view- 
ing angles close to 1/F as in our case; see Dermer 1995), in 
the comoving frame of a source at redshift z includes e.g. 
the CMB radiation field for large scale structures Uq MB — 
aT 4 (l + z) 4 r 2 (T = 2.7K is the temperature we observe 
now), and/or the photons of the Broad Line Region on sub- 
pc scales [/ BL r, etc. Thus the magnetic field required to give 
the observed L c /L s ratio is simply 

B = 5U ^jLzh^ G - (6) 

Using Eqs. 3 and 6 we can then constrain the values of F and 
B which minimize the total bulk power and at the same time 
ensure the required Compton over synchrotron luminosity 
ratio. Note that L B oc F 4 while L PiC oc S~ 4 . 

On the other hand, the Compton scattered radiation 
field can be dominated by synchrotron photons. Clearly, this 
would be favoured (for a given B) at low values of the bulk 
speed, such that the role of the external photons (i.e. their 
comoving energy density) is diminished in favour of the syn- 
chrotron ones (Sikora, Begelman & Rees 1994). In such sit- 
uation the Compton dominance is given by 

Lc _ 2 L s . . 

L s ~ cR 2 5 A B 2 ' () 
which leads to a value for the magnetic field 



■1/2 



B = 8.2 x 10" 



G. 



(8) 
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resulting in Lb oc F 2 <5 -4 and L p , oc T 2 . 

Note that the dependence of L PjC and L B with respect 
to the bulk Lorentz factor is radically different for these two 
possible origins of the high energy component. In the follow- 
ing we derive the minimum total power self-consistently, by 



taking into account which radiative process dominates for 
different values of B and V. 



3 APPLICATION TO PKS 0637-752 

Consider the specific case of PKS 0637-752, for which we 
have information on both the sub-pc and the extended jet. 

In Celotti et al. (2001) the spectrum of the X-ray knot 
WK7.8 was successfully fitted by adopting a one zone, ho- 
mogeneous synchrotron and inverse Compton model and as- 
suming relativistic bulk motion. As mentioned in the intro- 
duction, this was the key to explaining the strong X-ray 
emission (because the CMB is seen in the comoving frame 
boosted by the factor ~ T 2 and thus dominates over the lo- 
cally produced synchrotron energy density, in turn enhanc- 
ing the resulting X-ray emission). If the knot were not mov- 
ing relativistically, many more leptons would be required to 
produced the observed X-rays, since in this case: i) the flux 
would not be beamed and the intrinsic emitted luminosity 
would be larger; ii) the CMB would not be seen boosted and 
the inverse Compton process on external photons would be 
less favoured. 

Specifically, the model parameters adopted by Celotti 
et al. (2001) for the knot emission of this source assumed 
a bulk Lorentz factor F = 14, a viewing angle 8 = 5° and 
a source dimension R = 10 22 cm. Note that the emitting 
region is resolved in the X-ray images and thus R is not 
a free parameter. Also the viewing angle is limited by the 
VSOP observation of superluminal motion, but only to be 
8 < 6° at the pc scale; as radio maps do not show any 
bending of the jet between the pc and > 100 kpc scales, we 
assume the large scale jet to be at a similar viewing angle. 

The average {7), {-y 2 } derived by the spectral fitting are 
reasonably well constrained by the spectral slope and by in- 
terpreting the optical flux as due to synchrotron (this limits 
the maximum lepton energy 7 max ) and the X-rays as due 
to inverse Compton off the CMB. This latter constraint lim- 
its the minimum Lorentz factor of the distribution, 7 m i n , to 
values smaller than ~30, thus reducing the (usually large) 
uncertainty on the total power associated with this param- 
eter. Note that although the particle distribution is also a 
function of the assumed B and V, this dependence has an 
almost negligible effect on the power estimates In construct- 
ing Fig. 2 and Fig. 3, we have however taken into account 
the dependence of (7 2 ) on B and F by appropriately scal- 
ing in order to have the correct observed synchrotron peak 
frequency. 

From all the above parameters we estimated the pow- 
ers discussed in the previous section for the knot WK7.8 of 
PKS 0637-752. The main point to be stressed is that the 
total power (dominated by L p ) has a minimum for values 
of the Lorentz factor extremely close to that determined by 
the spectral fit fT = 11 vs F = 14). 

This applies of course for an ordinary e-p plasma as 
L p is calculated assuming an equal number of protons and 
emitting electrons. At this minimum the contribution of the 
Poynting flux is smaller (by a factor ~10) than L p and close 
to equipartition with the kinetic power associated with the 
emitting particles. Indeed in the case of an e* plasma, the 
total minimum is reached for F ~ 10 at the equipartition 
between Lb and L e . 
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Figure 2. Lower panel: The lepton and proton bulk kinetic pow- 
ers (L and L p respectively), the Poynting flux (Lb), the radiated 
[L r ~ (L' s + L' c ]r 2 ) and the total kinetic power (Ltot) as a func- 
tion of r for the parameters corresponding to the W7.8 knot of 
PKS 0637-752, as described in the text. The upper panels report 
the corresponding mass fluxes M p (in Mq yr — 1 ) and Dopplcr 
factors. The vertical dashed line corresponds to T = 15, the value 
used in the spectral fit. 



Given the relative freedom on the viewing angle and 
the fact that a decrease in 9 leads to a lowering of the to- 
tal required power, we tried to find independent estimates 
of Ltot- In particular we considered the correlation found 
by Rawlings & Saunders (1991) between the luminosity in 
the narrow line emission and the power required for jets to 
supply the extended lobe energy in their lifetime. From the 
luminosity in [OIII] (Tadhunter et al. 1993) we were able 
to estimate - following Rawlings & Saunders' prescription 
- the total luminosity in narrow lines I/nlr — 5 x 10 44 erg 
s _1 . According to the above correlation this corresponds to 
Ltot ~ 10 47 erg s _1 , although the correlation has (large!) un- 
certainties, and no energy other than that associated with 
the relativistic emitting particles in the radio lobes is as- 
sumed. 

Consequently we explored the possibility of a viewing 
angle smaller than 5°. In particular we report here the case 
of 9 — 3°, which roughly corresponds to a minimum power 
similar to that determined from the lobe energy content 
and line luminosity. The results are graphically displayed 
in Fig. 2. Similarly to the case already discussed the mini- 
mum value of the total power is reached for a bulk Lorentz 
factor extremely close to that inferred independently from 
the spectral fit (which required F = 15 and therefore almost 
the same B-field as before but a smaller intrinsic luminos- 
ity). In this case L p and Lb are closer to equipartition and 
about a factor 8 lower than for 9 = 5°. 

We also considered analogous estimates for the core 
(blazar) component of PKS 0637-752, assuming 9 = 3°, 
r = 18, R = 2.5 x 10 1B cm, Lblr = 10 45 erg s~\ Rblr ~ 
1 pc, B = 2.2 G, resulting in LT cx t ~ 0.1 and Ub ~ 0.2 erg 
cm -3 . Results are shown in Fig. 3. No consistent solution can 




r/3 

Figure 3. As in Fig. 2 for the parameters corresponding to the 
sub-pc core of PKS 0637-752, as described in the text. The ver- 
tical dashed line corresponds to T = 18, the value used in the 
spectral fit. 



be found at low F, where the dissipated power exceeds the 
total power carried by the jet. In addition, r ~ 18 is implied 
by the observed superluminal motion. Close to T ~ 18, the 
total power (once again dominated by L p ) presents a flat- 
tening (not really a minimum) at a value similar to that on 
extended scales and again for a V close to that inferred from 
the spectral fit (F = 18), and not far from equipartition of 
the magnetic field and lepton components. 

The two figures (2 and 3) also show in the upper pan- 
els the corresponding mass fluxes M p , estimated assuming 
one proton per lepton, and the Doppler factors. Again on 
both scales a minimum in M p is reached close to the mini- 
mum power condition and corresponds to values 0.1 — O.2M0 
yr . This can be compared with the accreting mass flux 
(note that at small scales entrainment is not likely to be 
substantial). If we consider a covering factor ~ 0.01 of the 
narrow line emitting gas, the corresponding ionizing (disc) 
emission requires an accretion flow of ~ IOMq^,^ yr -1 for 
a radiative efficiency of 77 = 0.1: only about 1 per cent of the 
mass would need to be channeled and accelerated in jets to 
ultra-relativistic speeds. 



4 DISCUSSION 

We have considered the implications of observed spectral 
quantities, namely the synchrotron (low energy) and inverse 
Compton (high energy) components which dominate the 
emission of extragalactic jets. Besides the total radiated lu- 
minosity, these quantities allow to estimate the correspond- 
ing bulk kinetic and magnetic powers as functions of the 
assumed bulk Lorentz factor. 

In particular it is possible to perform such estimates 
also for large scales jets which has been seen by Chandra. 
We have therefore analysed the case of PKS 0637-752, for 
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which broad band emission has been observed from both the 
blazar (sub-pc) and the extended jet. 

For the range of observing angles consistent with su- 
perluminal constraints we find that: a) the total power is 
minimized for relativistic bulk Lorentz factor larger than 10; 
b) the 'fitting' of the observed spectral energy distribution 
implies physical parameters which lie close to this minimum 
in the total power. These findings somehow lend support to 
the hypothesis of highly relativistic motion on the hundred- 
kiloparsec scale in jets, given the already high values of the 
estimated powers, which appear to be close to the limits of 
the proposed power production mechanisms in AGN. Note 
that the possible range of total power spans about three and 
the ratio L p /Lb about eight orders of magnitude. Further- 
more even the minimum power inferred (referring only to 
the non-thermal emitting jet plasma) is of order 10 47 erg 
s _1 , pushing (once again) the requirement for any powering 
mechanism. 

Within the acceptable range of viewing angles, a specific 
choice of 6 ~ 3° provides absolute values for the total power 
which well agree with both L to t derived on sub-pc scale and 
the jet power estimated from the lobe energy content (but 
note that this 'solution' implies to double the already large 
linear dimension of the jet). 

In other words it appears that jets are in the most ener- 
getically economic conditions for the given emitted radiation 
and spectrum, i.e. that there is a close link between the to- 
tal power, the dimension, the particle distribution and the 
intensity of the field. Despite of this the inferred radiative 
efficiency turns out to be only of the order ~ 3 x 10 -3 - 
10 -4 , for the nuclear and extended regions, respectively. It 
is worthwhile to point out that such low radiative efficiencies 
in the kpc-scale jet, and their higher values on smaller scales, 
may be a natural outcome of the internal-shock blazars sce- 
nario (e.g. Ghisellini 2000; Spada et al. 2001). 

In these conditions the mass flux in the jet corresponds 
to about one per cent of the mass flux accreting on the 
central source (for r\ ~ 0.1), as deduced from the luminosity 
of the narrow lines. 

It is relevant to comment on the possibility of an e 
dominated jet. Solutions close to the minimum total power 
for ultra-relativistic large (and small) scale jets can also be 
obtained. Furthermore such solutions appear to be intrigu- 
ingly close to equipartition values of L c and Lb- However, 
if the emission line luminosity is a measure of the lobe en- 
ergetics, then L c + Lb is too small to power the radio lobe, 
implying a contribution from the proton component, thus 
limiting the maximum pair content to 5-6 pairs per proton. 

On large scales the minimum power condition also cor- 
responds to a situation of quasi-equipartition between the 
magnetic and kinetic power transported in the jet. It is 
therefore interesting to ask whether there is a simple and 
robust way for the particle bulk motion and the magnetic 
field components to reach it. Theoretical studies of relativis- 
tic hydromagnetic flows (e.g. Li, Chiueh & Begelman 1992; 
Begelman & Li 1994 and references therein) have shown that 
rough equipartition between particle (kinetic) and Poynting 
flux is indeed a quite natural outcome for reasonable as- 
sumptions for the boundary conditions, with Poynting flux 
gradually accelerating matter. 

Clearly, various types of blazars present different spec- 
tral characteristics and radiated luminosities. We are already 



considering the implications for the total power involved ob- 
tained in the different cases (Ghisellini & Celotti 2001, in 
prep). A further analysis we intend to perform is on data 
relative to the extended emission from radio galaxies, in par- 
ticular to estimate the energetics which might be associated 
with a slower moving (external) layer. A final point to be 
stressed concerning the recent and future Chandra observa- 
tions is that although the procedure for the determination 
of powers is similar for the small- and large-scale jets, the 
large-scale powers are much better determined because the 
energy density of the external radiation field is well known. 



5 ACKNOWLEDGMENTS 

We thank the anonymous referee for useful comments which 
helped clarifying the paper. The Italian MURST is acknowl- 
edged for financial support (AC). 



REFERENCES 

Begelman M.C., Li, Z.-Y., 1994, ApJ, 426, 269 

Bridle A.H., Hough D.H., Lonsdale C.J., Burns J.O. & Laing 

R.A., 1994, AJ, 108, 766 
Burbidge G.R., 1959, ApJ, 129, 849 

Celotti A., Ghisellini G., Chiaberge M., 2001, MNRAS, 321, 
LI 

Chartas G., et al., 2000, ApJ, 542, 655 

Chiaberge M., Celotti A., Capetti A., Ghisellini G., 2000, 
A&A, 358, 104 

Dermer C. D., 1995, ApJ, 446, L63 

Fanaroff B.L., Riley J. M., 1974, MNRAS, 167, 31 

Ghisellini G., 2000, in Blazars Physics and Demographics, 
Urry M.C., Padovani P. eds ., Astronomical So ciety of 
the Pacific (ASP), in press ( |astro-ph/0011356| ) 

Gliozzi M., Bodo G., Ghisellini G., 1999, MNRAS, 303, L37 

Hardcastle M.J. Alexander P., Pooley G.G. & Riley J.M., 
1999, MNRAS, 304, 135 

Laing R., 1993, in Astroph. Jets, Burgarella D., Livio M., 
O'Dea C. eds., Cambridge University Press, Cam- 
bridge, 95 

Laing R.A., Parma P., de Ruiter H.R., Fanti R., 1999, MN- 
RAS, 306, 513 

Li Z.-Y., Chiueh T., Begelman M.C., 1992, ApJ, 394, 459 
Lovell .E.J., 2000, in Astroph. Phenomena Revealed by 

Space VLBI, Hirabayashi H., et al. eds., (Sagamihara: 

ISAS), 215 

Rawlings S.G., Saunders R.D.E., 1991, Nature, 349, 138 

Sikora M., Begelman M.C., Rees, M.J., 1994, ApJ, 421, 153 

Schwartz D.A., et al., 2000, ApJ, 540, L69 

Spada M., Ghisellin i G., Lazzati P., C elotti A., 2001, MN- 
RAS, in press ( |astro-ph/0103424| ) 

Tadhunter C.N., Morganti R., di Serego Alighieri S., Fosbury 
R.A.E., Danziger I.J., 1993, MNRAS, 263, 999 

Tavecchio F., Maraschi L., Sambruna R.M., Urry CM., ApJ, 
544, L23 

Tingay S.J., et al. 1998, ApJ, 497, 594 

Wardle J.F.C. & Aaron S.E., 1997, MNRAS, 286, 425 



© 0000 RAS, MNRAS 000, 000-000 



